A SIMPLE MODEL FOR NONEXPONENTIAL RELAXATION 

IN COMPLEX DYNAMICS 

A. Perez-Madrid 

Departament de Fisica Fonamental. 

Facultat de Fisica. Universitat de Barcelona. 

Diagonal 647, Barcelona. Spain 
The nonexponential relaxation ocurring in complex dynamics manifested 

in a wide variety of systems is analyzed through a simple model of diffusion 
in phase space. It is found that the inability of the system to find its equi- 
librium state in any time scale becomes apparent in an effective temperature 
field which leads to a hierarchy of relaxation times responsible for the slow 
relaxation phenomena. 
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Complex dynamics is the object of active research due to its implications in technology of 
materials and in several fields of scientific knowledge. At the physicochemical and biological 
level, complex dynamics is observed in glass-forming liquids, mechanical, dielectric and 
magnetic relaxation, amorphous semiconductors, pinned density-wave, protein dynamics, 
protein folding and population dynamics among others. The mechanisms underlying slow 
relaxation in complex dynamics still lack a clear and definitive elucidation. In the case of 
supercooled liquids and glasses, several experiments and computer simulations have been 
done which support the explanation of these relaxation phenomena in the framework of 
the energy landscape paradigm as the result of activated diffusion through a rough energy 
landscape of valleys and peaks [1]- [4]. 

To understand what these mechanisms are, we propose here a simple model to show 
a possible origin of nonexponential relaxation based on the idea of the energy landscape 
and nontrivial energy barriers. This model which consists of the diffusion in phase space, 



provides a direct link between the phase space dynamics and the slow relaxation of the 
functions of the configuration of the system in the corresponding energy landscape. The 
slowing down of the dynamics clearly appears as a consequence of the freezing of some 
degrees of freedom which takes the system out of equilibrium. This fact is indicated by the 
presence of an effective temperature field incorporating the information of the suppressed 
degrees of freedom and depending on the equilibrium temperature at the moment the quench 
was applied. 

Hence, we model the relaxation in the liquid as the Brownian motion of a test particle 
of unit mass in a potential. As it is well known, that physical situation is described by the 
Klein-Kramers equation which for simplicity we write in one dimension 



where p(x, u, t) is the probability density, x and u are the position and velocity of the test 
particle, 7 its friction coefficient and = k B T Q , with k B being the Boltzmann constant 
and T Q the bath temperature. Here, V(x) is a nonperiodic potential constituting a schematic 
representation of a rough energy landscape. The local equilibrium solution for the Kramers 
problem (1) is the Boltzmann distribution pi e q.( x i u ) ~ ex P { — P (1/2-u 2 + V(x))}. Thus, 
inherent to Eq. (1) is the existence of local equilibrium in phase space (x,u) [5], [6], and 
the approach to this state occurs at the bath temperature. Therefore, one concludes that 
at high temperatures the system possesses a large amount of energy to move through the 
whole phase space. 

We assume that a quench of the system freezes the traslational degrees of freedom taking 
the system away from equilibrium. Thus, this leads us to think of the partial decoupling of 
the probability density p(x, u, t) 



where the conditional probability <f> x (u, t; T Q ) describes a state of quasi-equilibrium [7], [8] in 
which the system, unable to equilibrate at the bath temperature, remains hanging. Assuming 
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that u is the fast variable, one concludes that the dynamical processes in the system are 
associated to configurational changes related to x which constitutes the slow variable whose 
probability density is p(x,t) = j pdu. Hence, by integration of Eq. (1) the time derivative 
of p is obtained 

d d r . . 

m p = -Txl updu (3) 

which defines the current J(x, t) = j updu. This current evolves according to 

which has been obtained from Eq. (1). Then, after partial integration of Eq. (4) and using 
the decoupling approximation (2), for times t ^> 7 _1 we obtain 

J(x,t) = -r^p(x,t)^-V(x) + k B ^-p(x,t)T(x,t,T )^ , (5) 

where r = 7" 1 and T(x,t;T ) = J u 2 (f) x (u,t;T )du is the second moment of the conditional 
distribution <p x (u,t;T ) which plays the role of an effective temperature which contains 
information on the frozen degrees of freedom. After defining the effective potential t) = 
V(x) + ksT(x, t;T Q ) Eq. (5) can be rewritten as 

J{x, t) = -D{x, t; T )faP(x, t) - rp(x, t)^(x, t) , (6) 

where D(x,t;T ) = Tk B T(x,t;T ) is the generalized diffusion coefficient. Thus, by substi- 
tuting Eq. (6) into Eq. (3) and taking r = 1 (i.e. rescaling the time t), this becomes the 
generalized diffusion equation 

l- p = ^^X)^p( X ,t) + pM^,t)) , (7) 

Note that the temperature field T(x,t;T ) introduces thermal barriers (or in others words 
non trivial activation energies) in the system [3]. The presence of these thermal barriers is 
a nonequilibrium effect that disappears when the system is in equilibrium. In addition, as a 
consequence of the elimination of degrees of freedom, the dynamics becomes non-Markovian 



and depends on the equilibrium temperature T a at the time of quench. By introducing the 
Massieu function 

e(», t) = io gP ( x , t) + j kBT( ^ ( , %) t)d x (8) 

the probability current Eq. (6) can be rewritten as 

d 

J(x,t) = -D(x,t;T )p Lqe Xx,t)—expQ(x,t) , (9) 

where pi. qe .{ x ,t) ~ exp|— / kB ^ x ^ -r^&(x, t)dx} is the local quasi-equilibrium probability 
density for which J(x,t) = or equivalently Q(x, t) = constant. 

If in a point £2 there is a sink where p(x2,t) = and a source at Xi so that a quasi- 
stationary current J(t) can be established in the system, thus, integration of Eq. (9) gives 
us the quasi-stationary probability density 

f x dv 

p(x,t) = J(t)p Lqe .(x,t) / — y . (10) 

Jx 2 D(y,t;T )p Lqe ,(y,t) 

On the other hand, a second integration leads to 

with n(t) = J^ 6 p(x, t)dx being the population at the left of the sink. Thus, we can formulate 
the following rate equation 

-n(t) = -K(t;x b )n(t) , (12) 



where 



/ f Xb f X2 dy \ 

K(f, x„) = d m , e .( X , t) I D(%t . To)pimM ) (13) 

is the rate constant, and Xb is the position of the barrier between x\ and x<i- The relaxation 
equation (12) admits the solution 

n(t) = n(t ) exp {- J* K(t'; x b )dt'} = n(t ) exp {-f(t)} (14) 



or 

n{t)=n{t ) exp{-/(t/r)} (15) 

which constitutes our main result, where n(t) = n(t/r) and the initial time t a is known in the 
literature as the waiting time. Here, the function f(t/r) might be interpreted as an algebraic 
function {t/r) 13 obtaining a stretched exponential behavior [9] or as a logarithmic function 
—a log At which leads to a power law behavior [10] that characterizes anomalous diffusion. 
In addition, another reading of our result is possible since assuming the form of f(t/r) 
is equivalent to assuming the form of the distribution of residence times ip(t) = — 4j.n(t) 
[11]. Thus, an in-depth analysis of the implications of the supposition of a nonexponential 
distribution of residence times usually performed in Continuous Time Random Walk models, 
reveals that this hypothesis might also be rooted in the nonequilibrium character of the 
dynamics inherent to the energy landscape picture as we have shown. 

Our previous discussion applies to fragile glasses. In the case of strong glassy systems, the 
potential is almost periodic. Thus, when the temperature is lowered, any configurational 
change is preceded by an equilibration in momentum space given by a local Maxwellian. 
Hence, the activation energies are independent from the temperature but with a certain 
random spatial distribution. Averaging over this distribution of energy barrier gives one the 
corresponding relaxation [4], [12], [13]. 

We conclude that despite its simplicity, the model we study contains the main features 
observed in the dynamical slowing down observed in a wide variety of complex systems. 
Hence, we think this work will contribute to the understanding of some aspects complex 
dynamics. 

This work was supported by DGICYT of the Spanish Government under Grant No. 
PB2002-01267. 
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